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We examine spinfoam cosmology by use of a simple graph adapted to homogeneous 
ff^ cosmological models. We calculate dynamics in the isotropic limit, and provide 

^^ the framework for the aniostropic case. The dynamical behaviour is calculating 

p^ transition amplitudes between holomorphic coherent states on a single node graph. 
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volume state, indicating that big bang type singularities are avoided within such 
models. 
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I. INTRODUCTION 

Cosmological models within the spinfoam framework serve a dual purpose [Ij: Their 
primary function is to form a proposal for extracting cosmological predictions from a full 
theory of quantum gravity. These models also perform a useful secondary role in forming a 
bridge between the canonical [2] and covariant [3] formulations of Loop Quantum Gravity 
(LQG). The covariant, or spinfoam, approach is a 'bottom up' construction - one predicates a 
quantum model and thence derives dynamics. As such the existence of a semi-classical limit 
and its agreement with the predictions of General Relativity are not a foregone conclusion, 
but rather must be examined within physical scenarios. This situation contrasts that of Loop 
Quantum Cosmology (LQG) [Ij, the application of the principles of LQG to cosmological 
mini-superspaces. 

Another important question that has to be clarified by spinfoam cosmology is whether 
physical predictions such as the resolution of cosmological singularities, a well trusted result 
in LQG [5] [6], can also be derived within this approach. The resolution of the big bang 
singularity, and its replacement with a deterministic bounce, is a key success of the canonical 
theory. It forms the basis of investigation of observable consequences of the theory [7j [8] . 
It is therefore a crucial test of the spinfoam approach that it reproduces these features in 
the ultraviolet sector. 

In |1] it was shown how to calculate the transition amplitude between two quantum states 
of gravity in the homogeneous and isotropic cosmological regime using a simple two-node 
graph (the dipole graph) at the first order in the vertex expansion. The main result of 
this work was to demonstrate that the used new spinfoam vertex amplitude (EPRL/KKL), 
P4T2]. together with some other ingredients, are adequate to derive a classical limit which 
can be identified with the Friedmann dynamics of an empty fiat, homogeneous and isotropic 
universe, i.e. Minkowski space. This result was further strengthened in [13], where a slight 
modification of the spinfoam vertex was utilized to implement a cosmological constant and 
derive a de Sitter universe as the classical limit. However, despite these original results 
being interesting, they exhibit a considerable deficiency, namely they fail to reproduce the 
curvature term 4"; which appears in the Friedmann equation. Such a term is expected, 
since the chosen graph is dual to a (degenerate) triangulation of the three sphere, (the 
closed topology, in which k = 1). The authors of [1] argue that this term might be recovered 
by taking higher orders of the spin approximation into account. However, this term appears 



in a more natural manner, as we will show in section [II G 

Our approach is the following: We will examine flat {k = 0) Friedmann-Lemaitre- 
Robertson- Walker (FLRW) models by use of a simple graph. This 'Daisy' graph consists of 
a single node which is both the source and target of three links. This graph can be thought 
of in two equivalent ways: In the first instance one has tessellated space by identical cubes, 
and so by symmetry opposite faces of a cube are identified, thus an outgoing edge dual to a 
given face is an incoming edge dual to the opposite face. The second instance is to consider 
the spatial slice to be a flat three-torus, with each link transcribing a compact direction. 
This is what separates this approach from the cubulation used in [T5] and also in [T5] . 

The second motivation is to provide the framework for investigating the more complicated 
case of Bianchi cosmologies. Since the inclusion of matter within the spinfoam paradigm 
has not yet been realized, one can only investigate FLRW models with trivial classical 
dynamics. In the anisotropic homogeneous systems comprising the Bianchi models there 
is a rich physical evolution even in the absence of any matter. These models have been 



examined extensively in LQC, both within the quantum framework [16], [T7], [18] and the 
semi-classical effective framework [19], [20], |2T]. 



II. RECAP OF THE THEORETICAL FRAMEWORK 

Let us briefly review the necessary theoretical input to make our ideas and calculations 
tractable, and fix our notation. Since we rely heavily on the theory as introduced in [1] we 
refer the reader to the original source or [3l [22] for a more detailed discussion. 

A. LQG and spinfoams 

The kinematical Hilbert space of LQG is defined as the direct sum of subspaces "Hr over 
all graphs F, embedded in a three dimensional manifold S. Since we want to work in a 
cosmological regime, describing just a finite number of degrees of freedom, it is sufficient 
for us to consider just one of these subspaces. This Hilbert space "Hr is defined on a graph 
r with L links and N nodes. Its elements are the spin network functions; gauge invariant, 
square integrable functions \l/ : SU{2)^ — )• C, {holonomy representation). Since gauge 
transformations act on the nodes N the Hilbert space "Hr is given by 

Hr = C\SU{2)''/SU{2f). (2.1) 

The name holonomy representation results from the circumstance that the SU{2) ele- 
ments hi are the holonomy of the Ashtekar-Barbero connection along the link /, i.e. 



hi = hi{A)=Vexpl A\ (2.2) 

with A = A^rMx". The components of A are given by A^ = F^ -|- jKl^ with F^ being 
the spin-connection, i^* the extrinsic curvature of S and 7 G M>o is the Barbero-Immirzi 
parameter. Thus the SU{2) elements hi contain the geometrical information of the quantum 
state ^(/i;)- 

Another representation, related to the former one via the Peter- Weyl transformation 
[21 [23], is the spin-intertwiner representation. In this representation the graph F carries 
spins ji G ^ at each link and invariant tensors in, called intertwiners, at each node. Those 
spins correspond to the spins of the unitary irreducible representations of SU{2) and the 
intertwiners belong to the S'[/(2)-invariant subspace /C„ = Invsuiijl'Hn], where "H^ is the 
tensor product of the representation spaces carried by the links meeting at the node n, 
"Hn = ®i(znl-lji. A general state in "Hr has the following structure in the spin-intertwiner 
representation 

^n.Shi) = f(g)^nj ■ (0V^''\hi)] , (2.3) 

where T>'^^'-\hi) is the 2ji + 1 dimensional Wigner matrix of the holonomy hi and the dot 
indicates contraction of indices. 

There are two interpretations of the role of the 3-dim. manifold S in which F lives: The 
original is that S is a spacelike slice at some coordinate time t. The spinfoam model would 
then define a transition amplitude from T-irc^t) ^o ^r(St+i)- This picture, making use of 



this time t, however, is not well suited if we want to preserve full 4-dini. diffeomorphism 
invariance. Thus, the second interpretation holds S to be a 3-dim. boundary of a 4-dini. 
spacetime region. The states in 'Hr(s) are thus not thought of as 'states at some time', but 
rather as boundary states, [5]. 

The dynamics of these quantum states can be defined via the spinfoam formalism. Think 
again of a boundary state \l/ G Tir with F C S. A spinfoam lives on a 2-complex made up of 
vertices, edges and faces. A 2-complex can be seen as a discretization of 4-dim. spacetime 
and heuristically may be thought of as resulting form a spin network evolving in time. A 
spinfoam model now assigns an amplitude to the state \1/ in the following way 

(^i^)=En^^(^) (2.4) 

where the sum is over different spin network histories and Wy{a) is called the vertex 



amplitude. We will present its precise structure in section |III| where we again will follow 
closely p. 

B. Coherent states 

Coherent states are an important tool for the examination of the classical limit of any 
quantum theory. In this section we will summarize a few definitions about the coherent 
states for LQG. In particular we will use the Livine-Speziale coherent intertwiners ^24j as 
well as the coherent states in the holomorphic representation [231 ES] later in this work. 

The Livine-Speziale coherent intertwiners make use of the Perelomov coherent states 



for SU{2) such that the intertwiner in in (2.3) is replaced by a coherent intertwiner. A 
Perelomov coherent state for SU{2) takes the highest weight state \j,j) G 'H^^\ which is a 
coherent state along e^, and rotates it with a Wigner matrix V^^\hfi) such that it is coherent 
along another axis n. The element /i^ G SU{2) corresponds to the 50(3) element Rft that 
rotates Cz into n. Thus we obtain the coherent state \j,n) = V^^\hfi)\j , j) . Consider a node 
n which joins E links e together. A coherent intertwiner at this node n is now given by the 
tensor product of the coherent states coming from each single link. The gauge invariance of 
these states is achieved via group integration. 

E 

$n(ne)=/ d(?(g)P(>)(^)be,^e) (2.5) 

JsUi2) g^i 

The holomorphic coherent states are characterized by an element Hi G SL{2, C) given at 
each link of the graph F. They are defined by 

-^mihi)= [ dgnllKtigsii)hig;'Hi), (2.6) 

where Kt is the analytic continuation of the SU(2) heat kernel to SL{2, C) and the group 
integration again ensures gauge invariance. The heat kernel is given by 

Kt{a, B)= J2 (2J + 1) e-"*^(^+^) TT{V^'\aB-^)) (2.7) 

jeNo/2 



with a e SU{2), B E SL{2, C) and a, t e M>o. The SL{2, C) label Hi now allows for two 
different decompositions ^23]. The first one is the polar decomposition 

and shows clearly that Hi determines a point in classical phase space on which the 
coherent state is peaked, hi G SU{2) is the holonomy of the Ashtekar connection A^ and 



El G su(2) is the flux of the densitized triad E"-. Thus, a coherent state with label (2.8) 
corresponds to a classical configuration {Al^,E!l). 

The second decomposition of Hi uses two SU{2) elements h^^ and h^' which, analogously 
to the SU{2) elements of the Perelomov coherent states, correspond to the transformation 
of Cz into fii and n^. Furthermore, a complex number zi is used whose real part is associated 
to the extrinsic curvature and its imaginary part is related to the area that is pierced by the 
link / Ea. 



,3 



Hi = h^^e-"'^hzl, (2.9) 

We denote the real and the imaginary part of zi as zi = q + ipi and a^ is the third Pauli 



matrix. The relation between the two decompositions becomes clear by writing (2.9) in the 
polar decomposition. One finds that [23] 

hi = hri,e-'''^'^h-_i,, (2.10) 



Ei= / E = 8nGhjpin[-—. (2.11) 

Where fi is the face dual to the link / with area Ai = SnGh'-f pi/ti. 



C. Classical preliminaries 

We are interested in the applicability of spinfoam cosmology to homogeneous models, 
both in the isotropic and anisotropic cases. In this section we will establish the holonomies 
and the fluxes for such models. We assume our spacetime to be of the form Ai = M x E, with 
S being a homogeneous 3-space. Under the additional assumption of isotropy the metric of 
A^ can be given by 

ds'^ = -dt^ + aitfd^'^ (2.12) 

with dVt^ = dr'^/{l — kr"^) + r'^dO'^ + r"^ siv? 9 d(j? and k G {0,±1}. The parameter k 
distinguishes three different spaces with constant curvature, where we are interested in the 
closed {k = 1) and the fiat {k = 0) case. The fiat and closed universes are special cases of 
the Bianchi I and IX universes respectively, in which all scale factors have been identified. If 



we consider a universe without matter but just a cosmological constant A, the metric (2.12) 
evolution obeys the Friedmann equation 



fc _ A 
a J ' a^ 3 



+ 7^=o- (2-13) 



In the case of vanishing cosmological constant the only possible solution is a static space- 
time a{t) = const., where for k = one recovers Minkowski space. If A 7^ one obtains 



for k = 0, and under the assumption that a{t) ^ 1 also for k = 1, the de Sitter solution, 
a{t) =exp{±^/A/3t). 

If we drop the restriction to isotropic models we obtain a Bianchi I universe in the flat 
case, which is described by the following line element 

ds'^ = -dt^ + ai{tfdx^ + a2{tfdy'^ + a^{tfdz^ . (2.14) 

Considering again a vacuum spacetime (with A = 0), the three directional scale factors 
oi, 02, as have to satisfy 

aid2Ci3 + 020.10.3 + 03O1O2 = . (2.15) 

This equation is solved by the so called Kasner universe and is given by aj(t) oc f^K The 
Kasner exponents have to fulfill the conditions Yli'^l — Yli'^i — 1- From those conditions 
one deduces that one exponent has to be negative, while the other two are positive which 
leads to a contraction in one direction and an expansion in the other two. 

Now, in order to specify the holonomy and the flux, we need the Ashtekar connection and 
the corresponding densitized triad. For that we will use the results provided in [HI I2SH2H] • 

In a general, i.e. non-cosmological, setting the Ashtekar connection is given by A^ = 
r^ + 7-?^*, with fC* beeing related to the extrinsic curvature 

Ki = e''Kab = y'C^a_)hab, (2.16) 

where the e^ are co-triads, such that the spatial metric can be expressed as hab = dije^e^. 
The connection coefficients T\ are calculated via contraction of the spin connection F^ = 
— ^e'^^'^Oajk, which is given by 

e:^ = -e^{d,el-rUe^ . (2.17) 

V^jj is the Levi-Civita connection compatible with hat, expressed in terms of the co- 
triads. However, using the framework of invariant connections on principal fibre bundles, as 



explained in [26j, simplifies the tedious calculation of A]^ via (2.16) and (2.17) enormously. 
Now, a Bianchi model is a symmetry reduced model of general relativity by a symmetry 
group S, which acts freely and transitively on S. If E is invariant under the action of S 
it is an homogeneous 3-space and a connection can be decomposed as A^ = (p^u^, with 
left invariant 1-forms wf and constant coefficients d>\. A further reduction, which leaves 



us for example with the three gauge invariant degrees of freedom in (2.15), is achieved by 
diagonalizing 0}. This has the effect that we can write 



A\ = c^k)A'k^: , (2.18) 

with A^ G SO (3), [S]. Using the left-invariant vector fields Xf, dual to the 1-forms u^, 
allows us to decompose also the densitized triad as 

i?f = pfx; = p^Af X^ , (2.19) 

where the second equality results again from a diagonalization of pf . These six coefficients 
{ck,P^), K = 1,2,3 now span the phase space of our reduced homogeneous model with the 
symplectic structure ^6] 

{c/,/} = ^75/- (2.20) 



If we expand the co-triads as ejj = e(x)A^a;f , with arbitrary ex G M, we get the following 
relations (no summation) 

P^ = l^iJK ejexl sgn(e/) . (2.21) 

With these simplifications the connection components F^ = r(/)A}(I;^ are given by, (no 
summation, even permutation of {1,2,3}) 

The n^ characterize our Bianchi model, we have for example n^ = for Bianchi I and 
n^ = 1 for Bianchi IX. Thus, we see that the Bianchi I models have vanishing spin connection 
r^. The extrinsic curvature is given by Kj = |e/, [5], where the dot indicates a derivative 
with respect to the coordinate time t. Now, we find the following results for the Ashtekar 
connection in the homogeneous setting 

A^ = c^j)A^,ui = (r(,) + |e(,)) A^jui . (2.23) 



In the isotropic case we have p^ = p"^ = p^ and (2.22) gives us Fj = in the flat case 
(Bianchi I), whereas we get F/ = | in the model with positive curvature (Bianchi IX). We 
can thus write c = ^{k + 'ye), k G {0,1}. For the anisotropic (Bianchi I) model we get 
ci = lei. 

Before we apply this formalism to our one-node graph let us make the following observa- 
tion. In [Ij it was shown that the holomorphic transition amplitude between two homoge- 
neous and isotropic quantum states, which are supposed to correspond to a curved geometry 
(fc = 1), is given by 

W{z)=NzeyiY>i-^ ■ (2-24) 

Following the reasoning in [13] the main contribution of W{z) is obtained when the real 
part of z^ vanishes and its imaginary part is proportional to vr/, / G Z. Now, if we use the 
correct relation between c and the metric variables, i.e. c = Re(2;) = \{k + •ja), instead of 
just c = d we can reproduce the correct Hamiltonian constraint. Therefore, we require that 
the real and the imaginary part (which doesn't contribute anyway, if we consider 1 1^(^)1) 
vanish. Thus, we get from z'^ = {c + ipY 



c^-p^ = 0. (2.25) 

However, the p^ term will disappear if we consider the proper normalized amplitude as 
done in [13] or [29]. Thus, we find c^ = and 



c' 


1 

~ 2 


/I -f^a^' 


\ 

= 
/ 




1 

= 2' 


l+7«) 


-i(- 


2 • 2\ 

7 « ) 




= c - 


-i(- 


7^d^) 






=^ 


-> 


- I'a') 


= 



(2.26) 
(2.27) 



Scaling of a and multiplication by a gives us 

-a^a + a = Q (2.28) 

which is, up to constants, the right Hamiltonian constraint for a curved FLRW universe. 



We have already mentioned the definition of the holonomy in (2.2). Now, let us define 
the flux. If we denote the link along which we evaluate the holonomy by / then Si denotes 
a surface pierced by /. One says S is dual to /. The flux of the electric field E = E°'T^Xa 
through a surface Si is given by 



E{S) = / {*Eyn^ , (2.29) 

Jsi 

where * denotes the Hodge dual, which converts our vector E into a 2-form, (dim(S) = 3), 
and n^ = nV* is a 5u(2) valued scalar smearing function [30]. We will use the definition 

(*E) = {*EyT^ = {*E)i^^^dx''' A dx^H^ (2.30) 

with 

i*E)i,^ = Baa.a.El (2.31) 

These definitions will become necessary especially for the anisotropic case, when we ex- 
plicitely have to calculate the Ashtekar connection and the flux for our model. 

III. OUR MODEL 

We want to calculate in this section the transition amplitude between two fiat, homo- 
geneous and isotropic universes. Therefore, we think of a cubical partition of 3-space S. 
The dual graph (with toroidal topology) is then given by the Daisy graph, see FIG.([I|. We 
will work in the one-vertex spinfoam expansion and we will use the coherent states in the 
holomorphic representation. The original motivation for this graph was its potential ap- 
plicability to anisotropic cosmological settings and therewith a physically more interesting 
situation, a problem we will tackle in a followup paper. In this first paper we will restrict 
our attention mostly to the isotropic case and see that this one vertex graph is sufficient to 
reproduce the original result of [1]. 

A. The Setting 

First, let us recall the definition for the vertex amplitude to specify the dynamics. Despite 
it being shown in [31] that there exist additional 2-complexes with the same boundary graph 
already at the one vertex level (for the dipole graph) we will consider just one spinfoam in 



the sum in (2.4). We consider the spinfoam that simply connects the two graph vertices 
with a single spinfoam vertex. The one vertex spinfoam expansion {v = 1) leads to the 
factorization of our amplitude Wv{a), [H [32] • With these simplifications and following [T] 
the transition amplitude between an initial and a final geometry is given by 

W{^f, ^i) = Wf{Hi{zf)) W,{Hi{z,)) , (3.1) 



where the single factors are 

W{z) 



QN 



dGnl[Kt{Hi{z),G,(^i)G 



t(/)) 



(3.2) 



where Kt{Hi{z),G) is defined in (2.7) and G is SO{A) for the Euchdean theory and 
SL(2, C) for the Lorentzian theory, respectively. As was explained in [33] for the Lorentzian 
case we will neglect one integration so that W{z) does not diverge. Furthermore, since our 
graph has just one node we find that source and target node of each link are the same, i.e. 
s{l) = t{l), thus leading to Gs{i)G~,i = I. 




FIG. 1: Cube and Daisy graph 

A clear advantage of using this graph is its simple application in the homogeneous case. 
It allows us to explicitly calculate the SL{2, C) elements for our coherent states and with 
that provides helpful insights also for more complicated structures. We begin by calculating 



the 5'L(2,C) elements Hi{z) using the decomposition ( |2.9[ ) 

Hi{z 



uie 



2 u 



r,-l 



(3.3) 



where ui and ui are elements of SU{2). We have three links, h, I2, h and six normal 
vectors rii = Cx and ni = ~ex,n2 = iy and n2 = —Cy and n^ = iz and n^ = — e^. The normal 
vectors ni and hi are obtained via a 5*0(3) transformation of e^ and the SU{2) elements ui 
and Ui are related to these 5*0(3) transformations, cf. appendix A. Now, we have to bring 
the three 5*L(2,C) elements in the following form 



HiU) 



3 2 3 

g-JaiVg-i/3Vg-Ja2V 



(3.4) 



This means that we have to find the angles ai, 02 and /3 . Given the SL{2, C) elements 
Hi{z) in this form we are then able to represent their Wigner matrices for all j if we recall 



that the angular momentum operators Jx, Jy, Jz are given hy Jx = ^ , Jy 
in the j = \ representation. Hence, we get 

V^'\Hi{z)) 



and Jz 






tU) 



-ia2J; 



U) 



One finds the following angles (cf. equation (A14), (A15), (A16) in the appendix) 

Hi(z) 



«! = ^2 



TT 

2 



H2{z) : 
Hsiz) : 






TT 

Z 



012 
02 



/3 






Tl — Z 



/3 
/3 



: TX 





(3.5) 

(3.6) 

(3.7) 
(3.8) 
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We can now calculate the transition amplitude W{z) = 

(3.9) 



[ dcflY, dj e-^"'''^^+''> Tr (v^^\Hi{z))G 



where we have defined dj = 2j + 1 and G = yl'P^-' '-^ \GsG^ ^)Y in the case of Euclidean 
gravity, (G e 50(4)), or G = Ftp{7i,i)(G'^G,-^)F in the Lorentzian case, (G e 5L(2,C))_. 
For detail cf. [3j. (Despite GG~^ = I, because s(/) = t{l) as mentioned earlier , we keep G 
for completeness.) Lets start by calculating the trace for / = 1 

j 
Tr(l?(^)(i7i(^))G)= J2 {j,m\V^'\H,{z))G\j,m) 

m=-j 

= E e-^^"^^''^U^£,{n-z){j,k\G\j,m), (3.10) 

m,,k=—j 

where we have inserted a unit operator. This leads to a separation of the geometrical 
information, stored in Hi{z), and the gauge invariant contribution, given by (j, fc| G |j, ttt,). 
For simplicity we will neglect the gauge contribution later on. For a detailed way of dealing 
with the trace cf. [22] • 

For Z = 2 we get analogously 

j 
Ti (v^^\H2iz))G) = J2 {j,m\V^^\H2iz))G\j,m) 

m=-j 

= E e-^--d^£,{n-z){j,k\G\j,m) (3.11) 

m,,k=—j 



Tt (v'^^\H,{z))g) = E {j,m\V^^\Hs{z))G\j,m) 



and for / = 3 



m=-j 

= E e-'^"^{j,k\G\j,m). (3.12) 

m,k=—j 

We will now use the large volume approximation, i.e. for Im(z) ^ 1 the term with 

m = j dominates. 

j 

Tr {v^^\Hs{z))g) ^ E ^"''' 0' ^1 ^ l^'^') • (3.13) 

k=-j 

Now, how do we treat the links / = 1 and I = 21 One can argue, that due to the highly 
symmetric setting we should also use m = j for those cases. If we do so we can make use of 
the following asymptotic relation 



d^l(p)=(-iy-"^J (M] ^x 

^"^^ ^ ' \l {j +m)\(j -my. 



X [cos(/3/2)]^+'" [sm{/3/2)y-"' . (3.14) 
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So lets start with 1 = 1. For m = j we get 



J2 e-'^'-""^^^ 4?(vr - z) {j, k\ G \j,j) . (3.15) 

k=-j 

By making use of 

€.(^-^) = (-l)™-^'rfS-^)W (3.16) 

and (3.14) we can analyse c?^'^ (tt — z) and get 



4J(vr-.)=(-l)^-^rf^L,)(^)=rfijL,)(^) 



X [cos(;2/2)]^-'^ [sin(2/2)]^'+'= (3.17) 

For the trigonometric functions we get for Im(2;) ^ 2 

cos(2/2) ^ ^e-^^/2 ^ sin(z/2) f^ ^e-^"/2^ (3.18) 

which gives us the following expression 

[cos(;2/2)]^"^ [sin(;2/2)]^'+' ^ 



2/ V2 

-J [ly+^e-'^K (3.19) 



We can use this to obtain 

2i 



Tr {v^^\Hi{z))G\ ^ (-\ e-'("+^)^x 



E -""^ (-ly^^'^^-^^y (— If^iyy a^|Gb,j). (3.20) 

If we take the e~*2-5 term inside the sum we get e~*2('^+-7) (—1)-'+^ V^^, which is equal to 

2i 



Tr(2)(^)(/7i(^))G)^Qy 



e-*"^x 



E (-i)^^""^^^""\/(— Ifri)! o-'^i^i^-'^-)- (3-21) 

k——j V 
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Now let us make a few comments about the last expression. First, notice that we can 
reproduce the factor e~^^^ , which appears also in the original work [1] and is crucial for 
the derivation of the transition amplitude. The second point is that we can proceed by 
approximating also the second part, namely the sum over k, as k = j, which simplifies the 
result and the important thing is, that by doing this approximation we are not doing worse 
then the projection onto m = j in the original work. 

The calculation for I = 2 is identical to the case 1 = 1. If we now use the following two 
relations 

_J =gln{l/4)j ^ (_l)J+fe,i+fc = (_^)i+fc_ (3,22) 

and apply the approximation k = j we get (neglecting factors like (—1)^-') 
Tr (v^^\H,iz))G) ^ e-^^^+i°(i)^- (j, j| G |j, j) , 

Tr (v^^\H2{z))g) ^ e-^^+'"(i)^- {j,j\ G \j,j) , (3.23) 

Tr (v^^\H,iz))G) ^ e-^^^ (j, j| G \j,j) . 

Now, let us compare these results with the calculations of the original paper [1]. There 
the authors used a projection onto the highest spin state m = j and got a factor exp{—izj) 
for all links. Our calculation gives the same result using a simpler graph. Furthermore, due 
to our explicit calculation we see in a more precise way where this factor comes from. Now, 
the term {j,j\ G \j,j) can be neglected in our case because s{l) = t{l) and thus 

(j,j|G|j,j) = (j,j|Ft2)(7.,.)(G,(,)G-(,ipr |j,j) 
= {j,j\Y^V^'^^'^\l)Y\j,j) 
= ((7J,J),J,J|I|(7J,J),J,J) = 1, (3.24) 

where we have used the projection property of the unitary y-map, cf. [I]. In [1] the 
authors replaced the whole gauge contribution ((j, j| G \j,j)) from all four links by a factor 
^ based on a calculation done in ^^. 

So far we have assumed that all our z-labels are the same for each link, which is a con- 
sequence of the isotropic configuration we are considering. If we investigate the anisotropic 
case these labels are going to be different z = zi. 

B. Results 



Inserting the results for the traces given by (3.23) and (3.24) into the transition amplitude 



(|3.9|) we get 

3 



Wiz) = I J2 (2j + 1) e-^'^'^^^+^^-''^ I . (3.25) 

Furthermore, we have assumed that Im(z) = p ^ ln(l/4) ^ —1.38, which is justified 
since p > 0, and thus all three links give the same contribution. 
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Now we can either apply a gaussian approximation as was done in [I] , we can investigate 
the amphtude numerically or we calculate (3.25) explicitly using the Cauchy product, all of 



which yield the same result. The gaussian approximation gives the result 



W{z) = (2jo + 1)= 



where Jq is given by 



Jo 



2th 



Im(2;) 
Ath 



3(2tft+iz)^ 
J Sth 



(3.26) 



(3.27) 



In order to get meaningful results we now have to normalize the amplitude. For this we 
use the following expression for the norm of a heat kernel coherent state given on a single 
link [29] 

4v^e*/^ 1 



^, 






(3.28) 



t3/2 sinh(p) 2 

where we have taken just the leading order term with n = 0, cf. [29]. The small g in 
the above formula corresponds to our 5*^(2, C) element Hi{z) and the heat kernel ti me t is 
related to our t via t = 2th. A detailed analysis reveals furthermore that the p in (3.28) 
corresponds our |. 



We will start with a numerical analysis of (3.25). Therefore, we plot A{z), which we 
define as the absolute value of W{z) divided by the norm to the third, because our graph 
has three links. 

Ai.)Jm- (3.29) 



« 



We set h = \ and truncate the sum (3.25) at j. 



150 where one has to make sure 



that jo < jmax, (3.27) holds. If we set the heatkernel time t = 1 we get FIG.(|2| and with a 



r(:) 




FIG. 2: Normalized amplitude A{z) for the Cube with jV 



150 and t = 1 



heatkernel time t = 0.1 we see that the peak becomes sharper FIG.(pj). 

Now, recall that Re(2;) corresponds to the extrinsic curvature of our model and thus we 
find Re(z) oc a = for all volumes Im(z) oc a^. Hence, we find that our universe is static. 
As we would expect. We don't have any matter or a cosmological constant, nor anisotropies. 

What is remarkable now is the decrease of the amplitude towards small Im(z). To see 



this more clearly insert (3.27) into (3.26) and calculate A{z) using (3.28). The result is 



A{c,p) 



sinh 



P\ _3£ _3cl 

- e 2 e 8*1 
2/ 



(3.30) 
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r(:>" 




FIG. 3: Normalized amplitude A{z) for the Cube with jmax = 150 and t = 0.1 



Plotting A{0,p) we get the shape in the p-direction which shows us the drop off for small 
scale factors FIG.Q. 



AO, p) 




FIG. 4: The transition amplitude as a function of scale factor, for a typical fixed c, here chosen 
to be zero. We find the amplitude is not supported on p = 0, indicating that there can be no 
transition to singularity. 



What does this tell us about the quantum dynamics of our model? Recall the definition 
of the transition amplitude between two quantum states of geometry \l/j and \1//, (3.1 ). One 



finds that the main contributions come from those configurations corresponding to classical 
geometries, namely d = and large scale factors a. The remarkable result is now, that the 
transition amplitude decreases for small scale factors and has zero support on a = 0. This 
is a statement about the occurrence of singularities in our model in that it tells us that a 
transition to a singularity is ruled out dynamically. 

The transition amplitude now gives us the possibilities for quantum fluctuations from one 
scale factor to another one. It doesn't give us the dynamical evolution of our universe. The 
classical notion of dynamics, i.e. d, is encoded in the phase space coordinates. In the same 
sense as a transition amplitude in QFT does not give us a temporal information of the time 
evolution of a certain process. 

Certainly this result has to be strengthened by future investigations, where it re- 
mains to show that one can circumvent the large spin approximation. In fact, one 
can calculate the traces in (3.9) in our model explicitly, thanks to the fact that our 

(j,fc|G|j,m) = (j,A;|rtp(7.,.)(G,(,)G7(,i^)y|j,m) = 

5km holds. This way one can in principle avoid 



graph has just one node and thus 
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the large spin approximation and indeed finds that the amphtude is still peaked on 
c = Re(z) = 0. However, the shape along the p = lm.{z) direction changes, a problem 
which probably has to be solved by the use of a different normalization. 

IV. DISCUSSION 

We showed in this paper that within the spinfoam cosmology approach the Daisy graph 
is sufficient to reproduce the vacuum Friedmann equation for a flat 3-space in the isotropic 
setting and thus may also be useful for the investigation of anisotropic models. The Daisy 
graph is perfectly suited to relax the restriction to isotropic models just by using three 
different holomorphic labels Zi at each link. This will be our setting in a following paper. 

Furthermore, we showed how one can reproduce the missing curvature term in the curved 
model, described by the Dipole graph, without using higher terms in the spin expansion. 
The right dependence of the real and imaginary part of the holomorphic labels z will also 
be important for the description of the anisotropic model. It has been suggested that the 
curvature term arises as a 'higher order' effect in the prior models. However, this derivation 
applies even in the fiat case, which would contradict the agreement with classical dynamics 
at large scale factor. 

The main result of this work, however, is the statement about the avoidance of singu- 
larities in our model. The dynamics of the Daisy graph show zero support for a transition 
amplitude from a finite scale factor to zero, therefore the singularity itself is not accessed by 
dynamics. Since the model under consideration does not include matter terms, there can be 
no direct comparison made with the bouncing models of LQC It remains to be seen at this 
stage whether this result is an artefact of approximations made, or is a deeper feature of 
the full covariant dynamics. A singularity resolution theorem or results analogous to those 
of [B] would be a strong achievement for the theory. In future investigations it has to be 
clarified, if this result survives if we loosen certain approximations we used in its derivation, 
and in the following paper we will investigate this in the anisotropic case. 

A serious open problem that remains to be investigated within the spinfoam cosmology 
approach is how to treat higher orders in the vertex expansion and of course the issue of 
gauge invariance for graphs with more than one node. A key motivation behind the use of 
the daisy graph is that gauge invariance becomes trivial, a feature due to the symmetries of 
the spacetime under consideration. 
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Appendix A: Details of the calculations 



In this section we present the calculation of the 5'L(2,C) elements Hi(z), (|3.3). As 



explained in section 2.1 the normal vectors rii and rii are obtained via a S0{3) transformation 
of iz and the SU{2) elements ui and ui are related to these 5*0(3) transformations. We start 



16 



with e-, ^-)■ Cr and e^ t— )■ —6^- 



I = 4"')e. 
,0. 



TT 

2 



cos(0) sin(0) 

1 
— sin(0) cos(0)^ 

l\ /O' 
10-0 
-1 0/ \1 



^0^ 



(Al) 








^-^^e. 



TT 

2 



/ cos(0) sin(0)' 

1 
\ — sin(0) cos(0) 



^0^ 



^0 -r 

1 

1 



^0^ 



(A2) 



We calculate the two corresponding SU{2) elements for the 5*0(3) rotation matrix R 
with the formula [35| [36] 



{12 + a^a'^Rrs) ^ e^^^o^ 
u = T -, — , / G SU{2) 

{2^1 + TiR) 



(A3) 



(x) 



We get for R^y> and Rl 



,(-^) 






R(-^0 ^ yC-^O = ^ 



Analogously one calculates e^ t-)- e^ and e^ ^-!■ — e^ with the resulting SU{2) elements 



(A4) 
(A5) 



i?l^) 



-^ u 



(2/) 



i?(-^) -^ 4-y^ = T 



I f I -i 



V2\-^ 1 



(A6) 



(AT) 



Finally, we have to calculate e^ i— )■ e^ and e^ h-> — e^ but it is clear that Ry is given by 



and thus the corresponding SU{2) elements is 



(A8) 



(A9) 
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Now we have to connect each two SU{2) elements with one hnk. We connect those vectors 
who are co-hnear. Furthermore, we need 



e— Y 



e-'2 
e'l 



The inverse matrices are given by 






-.^^^(l-l 



V2V 1 

1 /I i 



v^ V^ 1 



{u^-^)y' = Ti 



(AlO) 



(All) 



(A12) 



(A13) 



where we have chosen (A13) corresponding to the inverse elements of u*^^-* and m*-^^ resp. 



The reason is that formula (A3) is not valid for a rotation of n about the x— or y— axis. 



We get for the SL{2, C) elements 



HAz) 



Ml e '2 u 

lA-1 

2 U 1 



r.-l 



^We— V («(■ 



e-*t 
e^t 



-l' 



1 1 



g '2 — e 2 — e 2 — 6 2 



-isin(f) -cos (I) 
cos (l) isin (l 



-z I sm I - I a +COS I - ) o- 



and for Z = 2 and / = 3 we get analogously 



H2iz) = U2 e-'^"-2 u~^ = u^y) e"^"^ [u^- 



y) 



r 



M cos I -I (J -sm I - ) cr 



(A14) 



(A15) 



H.(z) 



-iz— ~-l 
M3 e 2 u^ 



e-*2 
e*t 



«(") e- 



{u^-r' 



(A16) 
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